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Abstract 

Striped Turing patterns and solitary band and disk structures are constructed 
using a three- variable multiscale model with cubic nonlinearity and global con- 
trol. The existence and stability conditions of regular structures are analysed 
using the equation of motion of curved boundaries between alternative states 
of the short-range component. The combined picture of transitions between 
striped and spotted patterns with changing level of global control is in quali- 
tative agreement with the results of the computer experiment by Middya and 
Luss [15]. 
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I. INTRODUCTION 



Stationary patterns of chemical activity (Turing patterns), have attracted wide attention 
both as a key to understanding complex and unusual phenomena in chemical and electro- 
chemical kinetics and as a possible basis of morphogenesis and biological diversification |l]-[3] . 
Similar patterns, emerging as a result of spontaneous symmetry breaking have been known 
and intensively studied long before that in fluid mechanics 

Reproducing stationary chemical patterns in the laboratory under controlled conditions 
turned out to be a rather elusive task, as most observations, such as the BZ reaction, 
produced propagating chemical waves rather than stationary structures; immobilization of 
chemical waves was achieved only due convective effects |J or imposed spatial gradients in 
the Couette reactor |7|||. In recent years, novel experimental techniques were applied to 
obtain a variety of complex chemical patterns in gel reactors pHl"2f and on catalytic surfaces 
T3[ . A parallel progress in computing and visualization is evident in 2d patterns obtained 



in recent computer experiments ||T4| . |T5| 



Most common models of large-amplitude dynamics, such as the FitzHugh - Nagumo 
equation, favor propagating waves, rather than stationary patterns. A basic mechanism 
generating stationary inhomogeneities is a combination of a short-range activator and a 
long-range inhibitor with a sufficiently fast dynamics [T(J. Formation of chemical patterns 
can be understood analytically when characteristic spatial and temporal scales associated 
with different reactants and other dynamic variables, like temperature or potential, are 
widely separated; this is typical to heterogeneous reacting systems, due to great disparities 
in diffusion rates. This scale separation allows to use rational approximation techniques, 
first of all, the method of matched asymptotic expansions, to build up and investigate 
dynamically complex non-equilibrium patterns |P7| , |18 . 



The basic element of a non-equilibrium pattern in a multiscale system is a propagating 
front, or kink, separating regions of prevailing alternative stationary states of the short- 
range subsystem, usually characterized by a low or high level of the activator. More complex 
patterns are constructed by allowing other constituents, varying on a longer spatial scale, 
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to modify the motion of kinks and to mediate their interaction, thereby leading to the 
formation of bound pairs, trains, etc. Generically, the kinks are "fasf , i.e. propagate with 
a speed of the same order of magnitude as a "unit" speed that can be conjectured from 
the characteristic length and time scales of the short-range activator. Wave patterns (with 
relaxation oscillations at any fixed point) are constructed as trains of fast kinks separated 
by distances far exceeding their width. 

Stationary or slowly evolving structures can be "assembled" in a similar way but using 
as basic elements "slow" kinks propagating with a speed much less than unity. The slowing- 
down can be either a consequence of a symmetry that is only weakly broken due to the 
influence of long-range fields, or a result of a dynamic adjustment to changing levels of a 
long-range variable. A Turing pattern can take the form of alternating domains separated 
by immobilized kinks. 

Dynamics of kinks is sensitive to the influence of global variables which are present in 



some of the most common laboratory set-ups. In experiments with catalytic wires [19,20 



the global control is due to heating input regulated to keep the average wire resistance 
(i.e. essentially the average temperature) constant. If local kinetics is bistable, the global 
control of this kind causes spontaneous symmetry breaking on a global scale [[TTJ due to 
immobilization of kinks. In electrochemical experiments |22|j2l| or glow discharge devices 



24J, the global controlling factor is the conservation of the total current. In both catalytic 



and electrochemical studies, oscillatory local kinetics, that would normally lead to a pattern 
of propagating waves (or coupled fast kinks) in a distributed system, generates standing 
waves under conditions of global control. Still more complex dynamics of pulses and kinks 
was obtained in numerical experiments [f25[| . 

A different kind of global control is connected with the influence of a rapidly mixed 
gas phase in catalytic reactions. Recent model computations |26[] indicated that a global 
gas-phase oscillation was the factor causing transition to standing waves in the experiments 
of the group of Ertl [I3|. Global interaction through the gas phase is usually detrimental 
to pattern formation, as it tends to synchronize the reaction on the entire surface. Recent 
computations by Middya and Luss |L5| demonstrated a sequence of transitions from striped 
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to spotted patterns and, eventually, to restoration of spatial homogeneity with an increased 
level of global control. 

The main purpose of this communication is to find analytic criteria for such transitions 
in a reaction-diffusion system with separated scales. Our description of Turing patterns 
is mainly based on representing them as collections of "solitary" band or disk structures, 
although we shall also consider regular striped patterns. Different objects bounded by kinks 
can be indeed viewed as solitary even at moderate distances, since interaction of kinks decays 
exponentially with separation. 

The motion of a slow kink is influenced by its geometric form even when its curvature is 
small. This distinguishes them from fast kinks, which are insensitive to the geometry unless 
their curvature radius becomes comparable to their width. Solitary structures bounded 
by slow kinks are stabilized by a delicate balance between the geometric factors and the 
influence of long-range variables, and their existence and stability conditions are strongly 
dependent on their geometric shape. 

In the following, we shall consider only slowly evolving structures that are built up of 
slow kinks, and have a characteristic scale far exceeding the width of the kink. The basic 
three-variable model with cubic nonlinearity is formulated in Section [TT], and the equation of 
motion of a bent kink is derived in Section |T|. We proceed with the construction of various 
stationary structures in Section |iy], and the analysis of their stability in Section [V]. The 
combined picture of the existence and stability regions of solitary band and disk solutions, 
as discussed in the last Section, is in qualitative agreement with the results of the computer 
experiment by Middya and Luss [O . 



II. BASIC EQUATIONS 

A. Reaction-diffusion model with separated scales 

The influence of global control variables on the formation of stationary non-equilibrium 
patterns can be understood qualitatively within the common context of pattern formation 
in reaction-diffusion systems with separated scales. The common recipe for a stationary 



4 



pattern |16|-|T^1 calls for a pair of equations for two reacting and diffusing species: a short- 



range "activator" u and a long-range "inhibitor" v. 

ju t = e 2 V 2 u + f(u,v), (1) 
v t = V 2 v + g{u,v). (2) 

The equations are written using the characteristic time and length scales of the long-range 
component as basic units. The small parameter e is the square root of the ratio of diffusivities 
of the two species. It is supposed that the nonlinear function f(u;v), with v considered as 
a parameter, has three zeroes within a certain range of v, two of them, u = M ± (f), being 
stable, and one unstable. Under these conditions, Eq.(JJ) has propagating kink solutions 
that separate regions with prevailing "lower" and "upper" states. Everywhere except in an 
0(e) vicinity of the kink, the dynamics of v is defined by a closed reaction-diffusion equation 

v t = V 2 v + h ± (v) (3) 

with a discontinuous source function given by alternative forms h ± (v) = gfw^f), t>) in 
alternative regions separated by a kink. It is clear that most of the information carried by 
the function g{u, v) is irrelevant. Its only essential property is reflected by the disposition of 
null isoclines, i.e loci of zeroes of the functions f(u,v) and g{u,v) in the plane (u,v), which 
can be either synclinal, as in Fig. |l]a, or anticlinal, as in Figs. |I|b,c. 

When the disposition is anticlinal, stable global stationary states satisfying h (v) = 
may either disappear altogether, as in Fig. [L]b, or become excitable, as in Fig. |l|c. An 
excitable state u = u (v) retreats when placed in contact with an alternative steady state 
of the the short-scale equation u = u T (v) at the same value of v. An advancing front of 
the alternative state is arrested when the level of the inhibitor adjusts to shifts in the front 
position, and relaxes to the value v = v s satisfying the stationarity condition for the kink: 

ru + (v s ) 

I{v s ) = / f(u,v s )du = 0. (4) 

An immobilized kink-antikink pair forms a solitary structure on the background of an ex- 
citable state. Such structures were constructed with the help of a multiscale model with 



piecewise-linear kinetics [p7| , and, recently, detected experimentally in gel reactors fL2 |. A 



regular array of alternating kinks and antikinks forms a Turing structure, which is the only 
stable stationary state in the "oscillatory" case corresponding to the disposition of null 
isoclines in Fig. |I|b. 

Stability of immobilized kink structures under conditions when the level of the inhibitor 
adjusts momentarily to shifts in the position of kinks can be supported by the following 
qualitative argument. We take note that, when the disposition of null isoclines is anticlinal, 
as in Fig. |l|b,c, v is depleted in the lower state u~ and produced in the upper state u + , 
while, on the contrary, the lower state advances at higher, and the upper, at lower values 
of v. Consider, for example, a finite region with u = u~ immersed in a continuum with 
u = u + . If it starts to shrink, the level of v would raise and, as a result, the kinks bounding 
the region would be immobilized again (a more precise analysis, taking also into account 
the local gradient of v and effects of the kink curvature is given in Section |V|). The same 
argument proves that immobilized fronts should be unstable when the disposition of null 
isoclines is synclinal, as in Fig. |l|a, so that the lower state advances at lower, and the upper, 
at higher values of v. 



B. Global control 

The above qualitative arguments can be extended to the case when v is a global controlling 
variable. The model with global interactions can be viewed as a limiting case of the model 
with separated scales, which is obtained when the diffusivity of the long-range variable tends 
to infinity, so that its diffusional length far exceeds the size of the system. Then Eq. (0) is 
replaced by the integral equation 

v t = (g(u,v)), (5) 

where (. . .) denotes averaging over the entire reactive surface. The respective form of eq. @ 
is 

Vt = ah~(v) + (1 - a)h + (v), (6) 
where a is the fraction of the surface occupied by the lower state. 



The action of either a long-range or a global inhibitor leads to the formation of stable 
inhomogeneous structures in the case of an anticlinal (but not a synclinal) disposition of 
null isoclines. The only distinction lies in the characteristic size of the inhomogeneities 
(either a wavelength of a regular pattern or a width of an excited domain). This size is, 
generally, determined by the diffusional length of the inhibitor, and increases to become 
comparable with the size of the system in the case of a global control. Under conditions of 
global control, the size of inhomogeneities is not limited, and the system evolves to minimize 
the length of the boundary between the domains with prevailing alternative states. Each 
disjoint domain tends to become circular as its circumference shrinks due to the line tension 
of a bounding kink. Further decrease of the length of the boundaries is caused either by the 
coalescence of domains, or by dissolution of smaller and growth of larger domains under the 
constraint of a conserved total area. This process is similar to spinodal decomposition in 
equilibrium systems described by the Cahn - Hilliard equation ||28|| . Its driving force is either 



the attractive interaction between like domains, which decays exponentially at a distance 
comparable with the diffusional length of the short-range variable, or the line tension, which 
causes the kinks to move with a velocity proportional to the ratio of the same diffusional 
length to the local curvature radius. As in equilibrium systems, the "ripening" process 
becomes exceedingly slow as the prevailing size of the extant domains increases. 

Separation into domains with prevailing alternative steady states in an anticlinal system 



is the basis of the "global regulator" model of spontaneous symmetry breaking |21|. This 
mechanism works, for example, in a catalytic wire heated by electric current when the 
controls are adjusted to keep the average temperature within the unstable region between 
the lower (extinct) and higher (ignited) states. In this case, a simple suitable form of the 
model equations is 

ju t = e 2 V 2 u + F(u) + U, U = B(u s -(u)), (7) 

which corresponds to an anticlinal disposition. The short-range variable u can be identified 
as temperature and the global variable U, as the overall heating rate. The disposition would 
not change qualitatively in a more precise model taking into account the dependence of local 
heat dissipation on local temperature. 



The situation is different when the global control is caused by the interaction with a 
rapidly mixed ambient fluid phase that changes its composition and/or temperature in 
response to the changes of the averaged state of the catalytic surface. The simplest model 
suitable to this case is 

ju t = e 2 V 2 u + F(u) + U, U -U = B{U - (u)), (8) 

where U has now the meaning of concentration or temperature in the ambient fluid. In this 
case, the disposition of null isoclines is synclinal, and no symmetry breaking occurs. If other 
factors induce formation of Turing patterns or excitation domains on the catalytic surface, 
the global interaction of this kind acts to modify the emerging patterns, and may suppress 
them altogether when it is sufficiently strong. 



C. Three- variable model with a cubic nonlinearity 

In the following, we shall consider a minimal model including three variables: a short- 
range activator u, a long-range inhibitor v, and a global controller of the synclinal type 
U that suppresses pattern formation. Analytical studies of non-equilibrium patterns often 
employ equations with a discontinuous piecewise linear source function for a short-range 
variable, as this form of the source function makes computations most simple. Such a model 
was used, in particular, to construct stationary structures as immobilized pairs or arrays of 



kinks ||27|| . In any physical situation discontinuous source functions can only appear when 
there is some hidden short-range variable that has been removed using a procedure of the 
same kind that has lead us from Eq.(|2]) to Eq.(0). In the long-range equation, however, the 
position of the discontinuity is not free but has a dynamics of its own, dependent on the 
dynamics of the short-range variable in the kink region. If this dynamics is ignored, and 
a discontinuous function is used in the basic short-scale equation, it can lead to undesired 
consequences, like the divergency of the propagation speed near the discontinuity. 

The simplest physically admissible source function in the equation of u should contain 
a cubic nonlinearity; the dependence on both v and U can be linear. We shall concentrate 
upon the case when the symmetry breaking between the alternative states of the short-range 



variable is weak, and causes the kinks to propagate with a speed comparable to that induced 
by a weak curvature with a radius far exceeding the diffusional length of the activator. This 
will allow us to take into account geometric factors while remaining within applicability 
limits of the multiscale perturbation scheme. The simplest suitable form of the model is 



With U < 0(1), eq.(|]) is just the Fisher - Kolmogorov equation supplemented by an 0(e) 
symmetry-breaking term. The latter is of the right order of magnitude to balance kinetic 
and geometric effects in the kink propagation. Both parameters n and (3 are supposed to 
be positive, which corresponds to the anticlinal and synclinal disposition of null isoclines, 
respectively, for v and U. Under these conditions, spontaneous pattern formation triggered 
by v is suppressed by the global variable which imitates the action of the environment in 
a mixed vessel. In the absence of global control ((3 = 0), both steady states u± ~ ±1 are 
weakly excitable at \v\ < //, while at \v\ > // only one state is excitable. 

The dynamics of the global variable is supposed to be very fast, and to insure stable 
relaxation to the solution of Eq. (pTTj) . The global balance condition ([11]) is written in such 
a way that, provided f3 = 0(1), the value of U is large, and relaxes to an 0(1) level only 
when the upper state u + ~ 1 is prevailing, so that the fraction of the surface occupied by 
the lower state is at most of 0(e). Therefore stable patterns generated by this model at 
moderate values of (3 can only take the form of a solitary band or disk with u < immersed 
in a continuum with u > 0. The only globally stable stationary state is, to the leading order, 
u — 1, v = n — u, U = Uq. At (3 < 7, another pair of solutions exists, one of which is 
unstable and the other one is metastable. 



e 2 V 2 w+ (1 -u 2 )u + e(U -v), 



(9) 




(10) 



(11) 



III. EQUATION OF MOTION FOR A CURVED KINK 



Problems with separated scales are generally solved using matched asymptotic expansions 
in the inner region (in the vicinity of the kink) and in the outer regions where one of the 
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alternative equilibria of the short-range variable is approached. We shall derive the equation 
of motion for a slowly propagating and moderately bent kink as a solvability condition that 
insures the existence of a stationary solution of Eq.([I]) in an aligned frame comoving with 
the kink. 

The coordinate frame aligned with a planar curve is defined in the following way. Let the 
curve be parametrized by a coordinate rj, and its position in the plane defined as x(?y). The 
derivatives of this function define the unit tangent vector t = x'(r])/\x'(rj)\ and the curvature 
vector dt/dr] = nn, where n is the normal to the curve, and k is the scalar curvature. By 
convention, the £ axis has its origin on the curve and is directed oppositely to n. The 
coordinate lines £ = const are obtained by shifting the curve along the normal by a constant 
increment; evidently, this shift causes the length to increase on convex, and to decrease on 
concave sections of the curve. The length element is defined as ds 2 = d^ 2 + (1 + K^) 2 drf. 
Accordingly, the Laplacian is expressed as 

v 2 = (i + ^y 1 [%(1 + + d v (l + k£)-%] . (12) 

The aligned system is well defined only sufficiently close to the curve, due to a singularity 
developing on the concave side at a distance of 0(n~ l ). If the curvature radius far exceeds 
the the diffusional range of the short-scale variable u, the aligned system remains regular in 
the region where u changes between the alternative equilibria u = u±. In the vicinity of a 
kink, the coordinate £ has to be rescaled by the factor e to accommodate the rapid change 
of u. Rescaling the tangential coordinate is not necessary as long as the curvature radius of 
the kink is of 0(1) on the extended scale, and far exceeds the kink width. Adopting this 
scaling, we rewrite the Laplacian Eq.(|T2"D as 



V 2 = e~ 2 <9 2 + e" 1 /^ + d 2 - k 2 ^ + 0(e). (13) 

The time derivative is transformed in the comoving frame as dt — > dt — Oe -1 ^. The 
additional term matches the curvature term by the order of magnitude, provided the propa- 
gation velocity C is scaled as C = e 2 7 _1 c where c = 0(1). In the inner region, the long-range 
variable is constant across the kink in the leading order. Assuming a linear dependence on 
some combination w of the long-range and global variables, as e.g. in Eq. (H), we set in 

ID 



Eq. (|l]) f(u,v) = fo(u) + ewfi(u) where fo(u) is odd. Expanding the short-range variable 
as u = «o(0 + £u^(£, rj) + . . . yields, in the zero order 

<(£) + /o(«o) = 0. (14) 

This equation is verified by a stationary kink (antikink) solution satisfying the asymptotic 
condition Wo(0 ~~ > u± &t £ ~~ * i 00 (f° r a kink) or at £ — > =Foo (for an antikink). In particular, 
for the cubic nonlinearity in Eq.®, the kink solution is 

Uo(0 = tanh(e/V2). (15) 

The first-order equation is 



u 



(i) 



+ ti(u )uU + *(e, i,) = 0, (16) 



which contains the inhomogeneity 

*(e,77) = (c + «K(0 + A(«o)w. (17) 



The propagation speed is determined by the solvability condition of Eq.(|T6|), which implies 
the orthogonality of the Goldstone eigenmode u' Q (£) of the homogeneous part, corresponding 
to the translational symmetry of a stationary kink or antikink, to the inhomogeneity rj): 

/OO 
u ' (0*ti,vH = o. (18) 
-OQ 

At w — 0, the well-known result c = — k (curvature driven motion [29[j) is evident. The kink 
propagates along the curvature vector in such a way that a convex region occupied by one of 
the alternative equilibria shrinks. The lowest-order general equation of motion for a curved 
kink is written as 

c{rj) = —k(t]) =F bw(r]), (19) 

where 

/oo / poo 

wo(0/iMOK// K(0M- (20) 
-oo / J — oo 



1 1 



In particular, for the model Eq.(|9]), u' (£) = 2 1 / 2 sech 2 (^/y / 2), fi(u) = 1, w = U — v, and 
b = 3/V2. The upper sign in Eq.([19|) applies to a kink with the lower state prevailing at 
£ < 0, and the lower, to an antikink with the reverse orientation. 

The equation of motion Eq.(|T9|) is closed by defining the instantaneous position of the 
kink in a suitable global coordinate frame and computing the local curvature. The outer 
equation Eq.(|2|) can be then solved by substituting the quasistationary values u = u in the 
regions separated by the kink, and using Eq.(|l9|) as a kinematic boundary condition at the 
kink. 



IV. REGULAR STRUCTURES 



A. Solitary straight kink 



The simplest structure in a system with separated scales is a steadily propagating solitary 
straight kink. Using in Eq.(|]) the quasistationary values of the short-range variable u = 
u ± (t>), which can be computed in the leading order in e using the roots u± of fo(u), and 
transforming to the moving frame we rewrite it as 

e 2 7 _1 c f x + v xx - v ± ix - v = 0, (21) 

where /1 = ^[g (u + ) — g (u~)], v = — |[(? (n + ) — g (u~)]; the same form is obtained directly 
when the model equation (^) is used. The small parameter in the first term appears when 
the propagation velocity is scaled as specified in the preceding Section. Provided 7 ^> e 2 , the 
propagative term is small, and will be further omitted. Then the leading-order unperturbed 
solution satisfying the continuity and smoothness conditions at the kink is 

u\e x — 11 — v at x < 0, 

(22) 

n[\ — e x ] — v at x > 0. 

When the short-scale equation @ applies, the stationary propagation speed, defined by 
Eq. ( |19|) is c = —b[U — v(0)} = —b(U + v). The kink is at rest when the fraction of the 
surface occupied by the lower state, a, is a = (e/2/3)(Z7 + v)- Since this ratio is small at 
(3 = 0(1), it is likely that boundary effects have to be taken into account in this case. 
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B. Striped pattern 



A regular striped pattern in the infinite plane with the period 2L = 2(a + + a_) is 
formed by alternating straight-line kinks at x = 2nL and antikinks at x = 2(nL + a + ). The 
leading-order stationary solution for the long-range variable is 



V = yU 

V — /i 



smha_|_ 

2 cosh (x — 2nL + a) — 1 

smh L 

sinh a_ 



1 — 2 cosh (x — 2nL — a + 



sinh L 



v at 2(nL — a_) < x < 2raL, 
z/ at 2nL < x < 2( , nL + a + ). 



(23) 



The equilibrium condition, that fixes the relation between the lengths a±, is v(0) = v(2a + ) 
U. A short computation yields 



U + v — /isinh(a + — a_)/ sinh(a + + a_). 



(24) 



Using also Eq. (p^lf ) we obtain 



U + v = n 



sinh(a + — a_ 



2/3a_ 



sinh(a + + a_ 
Evidently, a + ^> a_, unless /3 < 0(e). 



e(a_i_ + a_ 



(25) 



C. Solitary band 

A solution in the form of a solitary band with u < immersed in a continuum with 
u > can be obtained from Eq. ( PB"[ ) in the limit a + — > oo. The solution, that satisfies the 
continuity and smoothness conditions at the kink and the antikink located at x — ±a, is 



fj,(2e a coshx — 1) — v at |se| < a, 
fj,(l — 2e~\ x \ sinh a) — v at \x\ > a. 

Some typical profiles of v are shown in Fig. ^|a. The stationarity condition is 



(26) 



v + U = fie 



-2a 



(27) 



Solitary structures naturally appear at (3 — 0(1). Since the interaction between distant 
kinks is exponentially small, one can also envisage a pattern consisting of N well separated 
parallel bands. When such a structure is placed on a plate with the width L/e, the fraction 
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of the surface occupied by the lower state is a = 2eNa/L. Then U = Uq — 2/3ia, where 
(3\ = 2Nf3a/L. It is easy to see that Uq can be compensated by a shift of u; further on, we 
shall set it therefore to zero. Then the stationarity condition Eq. (p?D can be written in the 
form 

v - 2/3 1 a = /ie" 2a , (28) 

If this equation is solved graphically, stable solutions can be obtained only when the absolute 
value of slope of the straight line representing the left-hand side is less than that of the 
right-hand side at the intersection point. The necessary condition for the existence of such 
solutions is fti < \i. Within the parametric range /i > (3\ > 0, \i > v > v c = /3i[l — ln(/3i//i)], 
there are two solutions, of which one corresponding to a smaller width is stable according 
to the above criterion. The largest possible width, achieved at v = u c , is a = |ln(/z//3i). 

D. Solitary disk 

In 2d, the kink should acquire a circular shape due to the "line tension", and the sta- 
tionary solution depends only on the long-range radial coordinate r. The stationary outer 
equation is 

v rr + r~ x v r — v — v ± fi = 0, (29) 

The stationary solution for a disk with u < immersed in a continuum with u > 0, satisfying 
the continuity and smoothness conditions at the kink located at r = a, is 

{fi[2aK 1 (a)I (r) — 1] — v at r < a, 
(30) 
fj,[l — 2aIi(a)K (r)] — v at r > a, 

where I n , K n are modified Bessel functions. The identity Ki(a)I (a) + K (a)Ii(a) = a" 1 
should be used when the boundary conditions are checked. Some profiles of v are shown 
m Fig. |b. One can notice that the value of the lonEf-range variable at the kink location is 
shifted upwards to compensate the curvature-driven shrinking action; this shift grows with 
the decreasing disk radius. 



14 



Using v(a) and k = a 1 in Eq.(|19D yields the stationarity condition 

U + u= -{ba)- 1 + fw, [K^a) I (a) - K (a) Ji(a)] . (31) 

If a pattern consisting of N well separated disks is formed on a plate with the surface area 
(L/e) 2 , the fraction of the surface occupied by the lower state is a = enN(a/ L) 2 . Then 
U = Uo — 2tt (3N (a / L) 2 , (where Uo can be again set to zero), and the stationarity condition 
Eq. ( |3~ll) can be written, using the parameter (3^ = neNf3/L 2 , in the form 

va + b' 1 = 2(3 2 a 3 + /ia 2 [K t (a) J (a) - K (a) h(a)) . (32) 

If this equation is solved graphically (Fig. |3|), the slope of the left-hand side of this equation 
at the intersection point corresponding to a stable solution should be algebraically larger 
than the slope of the right-hand side. When it is so, the effective repulsive action due to 
the long-range field decreases with growing radius faster than the shrinking action due to 
the line tension. Whenever Eq. (pT|) has a single solution, this solution should be unstable. 
When an additional pair of solutions bifurcates at a tangent point, one of these solutions 
is stable according to the above criterion. Further stability analysis is found in the next 
Section. 

V. STABILITY OF REGULAR STRUCTURES 

The "naive" stability criteria in the preceding Section refer only to the stability to per- 
turbations that do not reduce the symmetry of the problem. In this Section, we shall 
investigate stability to asymmetric perturbations. Perturbations of this kind do not change 
the fraction of the surface occupied by either state, and therefore do not involve the global 
control variable. Throughout this Section, we shall assume that 7 is not exceedingly small, 
i.e. 7 > 0(e). Under this condition, the distribution of the long-range variable follows slow 
shifts of the kink position quasistationarily, unless at distances that are large even on the 
extended scale. 
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A. Stability of a solitary straight kink 



We shall now apply the equation of motion Eq. ([l9D to the problem of stability of a 
solitary straight kink to infinitesimal perturbations. The suitable global coordinate frame 
for this problem is a Cartesian frame (x, y) comoving with the unperturbed kink. Denoting 
the instantaneous position of the kink relative to its unperturbed position as x = C(y,t), 
we take note that, as long as the amplitude of the perturbation is much smaller than its 
wavelength, the normal vector, defining the direction of the kink propagation, is almost 
parallel to the x axis, and the propagation speed is expressed as ( t — e 2 7 _1 c. In the same 
approximation, the curvature is given by k = —( yy - Expanding also the long-range variable 
in the vicinity of the unperturbed position, we reduce Eq.(|19|) to the form 

(r = Cyy + b[V(0,y,r) + v' s (0)Q. (33) 

Here v' s (0) = fi is the stationary flux at the kink following from Eq. fl22|) , V is the perturbation 
of the long-range variable, and r = e 2r y~ l t is the rescaled time variable. 

The correction to v due to a weak perturbation of the kink can be computed in a simple 
way by observing that any shift of the kink position by an infinitesimal increment ( is 
equivalent to switching the sign of the source term in Eq.(|l(]) in a narrow region near x = 0. 
This approach is akin to the "singular limit eigenvalue" method applied earlier to the stability 



analysis of kinks in one dimension |30| . Although the perturbation £ is habitually presumed 
to be arbitrarily small for the purpose of linear stability analysis, the method is actually 
applicable also to finite perturbations of the kink position, provided they are small on the 
extended scale, i.e. ( < 0(e). 

At 7 > 0(e), the quasistationary form of the perturbation equation is applicable: 

V 2 V-V = 2rt{y,t)6{x). (34) 

In a standard way, we are looking for a solution in a spectral form 

C(y, r) = X (k, X)e Xr+iky , V(x, y, r) = ^(x; k, X)e Xr+iky . (35) 

The dispersion relation X(k 2 ) determines the stability to infinitesimal perturbations. Using 



Eq. (|35D in Eq.fl34|) yields 



?pxx - q 2 ip = 2{J>((y, t)6(x), (36) 
where q 2 = 1 + k 2 . The solution of this equation is 

^(x) = -/i X (A;,A)g- 1 e-^l (37) 

Using Eqs.(|22|), fl35|), ( |3?| ) in Eq.(p3|) yields the dispersion relation 

\ = -k 2 + bfi(l- q- 1 ). (38) 

The eigenvalues \{k 2 ) are always real, and, as expected, A vanishes at k = 0, which 
reflects the translational symmetry of the kink. The loss of stability to long-scale pertur- 
bations occurs when the derivative d\/d(k 2 ) at k — becomes positive. This happens at 
bfi > 2. Since the function A(fc 2 ) is convex, the loss of stability always occurs with growing 
\i in the long-scale mode. 



B. Stability of a solitary band 

The equations of motion for infinitesimal displacements C^, of the kink at x = a and 
the antikink at x — —a that confine a solitary band (Section [IV C| ) are written analogous 
to Eq. (|3~3"D . We take note that opposite signs appear in the equations of motion Eq. (pj|) 
for the kink and the antikink. Using the expression for the stationary fluxes following from 
Eq.(pBl), we obtain 

(39) 

It is advantageous to consider the symmetric and antisymmetric combinations of the dis- 
placements, C ± = |(C (1) ±C (2) ), that may lead, respectively, to zigzag or varicose instabilities. 
It is clear from Eq.(p9|) that £ are coupled, respectively, to the antisymmetric and sym- 
metric parts V T of the perturbation field V: 

(40) 

The quasistationary equations for both perturbation fields, V T , are identical: 







(?) 



Si 



'/// 



+ b 



±V(±a,y,r)+fi(l-e~ 2a )C 



C ± = c 



± 
w 



V^(a 1 y,T) + f i(l-e~ 2a )C 
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VV T - V T = 2^( ± (y, t)6(x - a). 



(41) 



Setting as before V ± = ^ e XT+iky , C ± = x ± e Xr+lky yields 

— 2[j,x + q~ 1 e~ qa coshgx at \x\ < a, 
— 2fix + Q~ 1 e~ q ^ coshga at \x\ > a, 
—2nx~q~~ 1 e~ qa sinhqx at |x| < a, 

— 2 sign(x) fXX~Q~ le ~ sinh qa at \x\ > a. 
The zigzag and varicose branches of the dispersion relation are given by 



ip (x) 



(42) 



(43) 



A + = -k 2 + 6/i 1 - e~ 2a - q- 1 (l - e- 2qa 

\- = -k 2 + fe/x [l - e~ 2a - q- 1 (l + e- 2qa )] . (44) 

Vanishing A - at k = reflects the translational symmetry of the band. Due to convexity 
of the function X~(k 2 ), the loss of stability first occurs at k = 0. The solution is unstable 
when the derivative d\~ jdk at k = is positive. The critical point of the long-scale zigzag 
instability is 

2/b 



(45) 



[l + 2a)e- 2qa - 1' 

The varicose branch, X + (k 2 ) lies persistently lower than the zigzag one. Thus, the stability 
limit is always given by Eq. (flop . The most dangerous perturbation corresponds to the long- 
wave limit k — > 0, i.e. q — > 1. The limit of the zigzag instability is shifted to higher values 
of n with a decreasing band thickness. 



C. Stability of a solitary disk 

In order to check stability of the stationary solution that has the form of a solitary disk 
with the radius a, as in Section [IV D| , we use the cylindrical coordinates (r, 0), and set r = 
a[l+(((f), t)], where £(</>, t) <^ 1. The curvature of the kink is expressed as k = a _1 (l— (— C/></>)- 
Expanding, as before, the long-range variable in the vicinity of the unperturbed position 
brings Eq.(|l9D to the form 



+ () + b a-V(O,0,r) + <(O)C 



(46) 
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where v' s (0) = 2fiaIi(a)Ki(a) is the stationary flux at the kink following from Eq. (|30|) . The 
quasistationary equation for the perturbation field V is 

VV-y = 2fiaC((f),t)5(r-a). (47) 
We are looking for a solution in the form 

C(0, r) = X ne XnT+m t V(r, 0, r) = ^(r; n)e x - T+in \ (48) 

where n is an integer. The equation of ip is 

1 / n 2 \ 

tfj rr + -lfj r - 1 + — = 2/iaXn5(^ ~ «)• (49) 



The solution of this equation is 

!—2ua 2 K„(a)I n (r) at r < a, 
(50) 
—2fia 2 I n (a)K n (r) at r > a. 

Using Eqs.(|30|), fl4"8|), ( |50| ) in Eq.(p3|) yields the dispersion relation 

A n = -{n 2 - l)/a 2 + bna[h{a)K x (a) - I n {a)K n {a)}. (51) 

As expected, A n vanishes at n = 1, which corresponds to a shift of the disk without 
deformation. One can check numerically that, of all angular modes, the quadrupole mode 
n = 2, deforming the disk into an ellipse, is the most dangerous one. At large a, this can be 
seen by using the asymptotics of Bessel functions to compute the asymptotic value of /i at 
the instability threshold 

(n 2 - l)/b 



fi c {n) 



2 



_9_ 

61 



(52) 



i[r(| + n)/r(-| + 

where T(x) is the gamma-function. The asymptotic limiting value increases with n at n > 2. 

D. Stability of a striped pattern 

A quasistationary perturbation V of the long-range field in alternating regions of a striped 
pattern with the period 2L (Section |IV B|) due to displacements of respective kinks at x = 



2nL by Q 1,n ) and antikinks at x = 2{nL + a + ) by C( 2,n ) obeys the equation 
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V 2 ^ - V = 2/i^(-l)" \C ihn) 5(x - 2nL) - C (2 ' n) 5(x - 2nL - 2a, 



(53) 



The equation of motion of the kinks is obtained, analogous to Eq.(^), by using in Eq.(^9|) 
v = v' s (2nL)(( n ^ + V(2nL,y,r), where v' s (2nL) = J = 2/zsinha + sinha_/sinhL is the 
stationary flux following from Eq.(^); at the antikink locations, x = 2(nL + a+) , v' s (2nL) = 
—J, and the respective values enter the equation of motion with the reverse sign: 



C^' n) = CS' n) + b \ JC (1 ' n) + V(2nL, y, r) 



(2,n) + b 

yy ^ 



j C (2,n) - V (2nL + 2a + ,y,r) 



(54) 



Presenting the perturbation V in each successive stripe with a prevailing positive or negative 
state of the short-scale variable as V^ ,n ' = il)^' n )e XT+tky , and kink displacements as £^> n ) = 



■yb » e ^r+iky ^ j _ 2 ; we write the equation of the spectral components 



(55) 



The matching conditions, taking into account localized perturbations due to displacements 
of respective kinks and antikinks, are 



^(l,n) = ^(2,n) ) ^(l,n) _ ^(2,n) = 2; u X ( 1 >™) at X = 2nL, 

^(i,n) = ^(2,n+i) ) ^,(i,n) _ ^,(2,n+i) = 2fi X (2 ' n) at x = 2(nL + o+). 
The solution of Eq.([5^) is expressed as 



(56) 



(57) 



We shall define the coefficients in the above expression, as well as the displacements at 
particular locations, through respective generating functions: 



aj{z) 



]T Ab> n ) e {iz+2qL)n , pj(z) = ]T B^ n) e {iz ~ 2qL)n , 



7i(*)= E X 0>) e izn , J = 1,2. 

n=— oo 



(58) 

Multiplying the matching conditions Eq.(|5B]) by e lzn and summing over n yields a set of 
equations for Oj(z), (3j(z): 



2n 



e 



2qa 



ai(z) - a 2 (z) + (3i(z) - (3 2 (z) = 0, 
«iO) - a 2{z) - (3i(z) + (3 2 (z) = 2[iq~ 1 j 1 (z), 
atiz) - e- {lz+2qa - ] a 2 {z) + e^+P^z) - e - (iz ~ 2qa - ) f3 2 (z) = 0, 



? qa +a x (z) - e- (l2+2ga - ) a 2 (^) - e~ 2qa+ p 1 (z) + e^- 2 ^ f3 2 {z) = 2 m ~ x l2 {z). (59) 



Two more equations are obtained in the same way from Eq. (|54|) : 



(\ + k 2 -bJ) ll (z) = a 1 (z)+f3 1 (z), 

(X + k 2 - bJ) l2 {z) = -e 2qa + ai (z) - e- 2qa +p 1 (z). (60) 

The general dispersion relation, obtained as the condition that Eqs. (|59|), (|60|) have 
a non-trivial solution, is quite cumbersome. As in Section |VB| , it is convenient to work 
with zigzag and varicose modes, that correspond, respectively, to like and opposite signs 
of displacements of adjacent kinks and antikinks. These modes, obtained by introducing 
the symmetric and antisymmetric combinations ^ = Mji ± 72), decouple, however, only 
at z = 0. The symmetric (zigzag) and antisymmetric (varicose) branches of the dispersion 
relation at z = are computed as 

' sinh a + sinh a_ sinh qa + sinh ga_ \ 



A + = -hr + 2bfi 



sinh L q sinh qL 



2 ; ( sinh a + sinh a_ cosh qa + cosh ga_ 



A' = -k 2 + 2b» _ -7- — . (6i) 

y sinh L g smh gL y 

The varicose branch A _ (A; 2 ) lies persistently lower than A + (/c 2 ), and therefore can be 
disregarded. Vanishing A + at k = reflects the translational symmetry of the entire pattern. 
The loss of stability in the long-scale zigzag mode occurs when the derivative d\ + /d{k 2 ) at 
k — vanishes. 

A simpler form of the dispersion relation is obtained at v = when the striped pattern 
is symmetric, with a + = a~ = L. The zigzag branch of the dispersion relation is 

q^ 1 tanh qL 



A = —k + bfj, 



tanh L 



(62) 



1 - sin 2 f (l - tanh 2 qL) 

One can clearly see here that A decreases monotonically with growing z, so that when /1 
increases the loss of stability always occurs in a concerted long-scale zigzag mode z = 0. 
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Again, the bracketed expression in Eq.(|62]) is a convex function of k 2 , and therefore the 
loss of stability should first occur at k — 0. The instability limit shifts to larger /i as the 
wavelength decreases, as shown in Fig. [|. In the other limit, a + ^> a~ , Eq. reduces to 
the dispersion relation for a solitary band, Eq. 

VI. DISCUSSION 

Although both the zigzag instability of the solitary band solution and the quadrupole 
instability of the solitary disk are not coupled to the global control parameter /3, the stability 
boundary in the parametric plane (//, v) depends on the global parameter as well because of 
the parametric dependence of the stationary solution. The existence and stability regions of 
the solitary band solution are shown in Fig. [|. In the absence of global control, the solution 
exists in the entire triangle fi > v > 0, but suffers zigzag instability to the right of the 
dashed line departing from the // axis. Both the existence and stability limits shift upwards 
as the parameter /3\ increases. 

The existence and stability regions of the solitary disk solution are shown in Fig. |^. At 
(3 2 = 0, the tip of the existence region penetrates into the half-plane v < 0; at // > the 
lower boundary of the existence region coincides with the axis v — 0. Two other cusped 
curves, drawn for j3 2 = 10~ 2 (Fig. ||a) and (3 2 — (Fig. ||b), bear witness that the 

existence conditions are very sensitive to small changes of the global parameter (3 2 . This 
sensitivity is, indeed, essential, since the definition of the parameter /3 2 includes the small 
parameter e, and therefore (3 2 (or, at least, the ratio f3 2 /f3{) should be small. The limits of 
the quadrupole instability are shown in the same picture by dashed lines (at (3 2 = 10 -3 , the 
instability limit is outside the depicted region). 

The regions where either or both band and disk solutions can appear are shown together 
in Fig. 0. The coexistence region is most extensive when there is no global control (Fig. |7|a,b). 
One can recall that coexistent striped and dotted patterns (which can be viewed as arrays of 
band or disk structures) were also found experimentally Jll| . A narrower coexistence region 
is seen in Fig. [7]c, where the region of stable solitary disk solutions at f3 2 = 10~ 3 , bounded 
by the solid cusped line, is shown together with dashed lines representing the lower bounds 

22 



of the existence of stable solitary band solutions at several values of (3\. The coexistence 
region shrinks with increasing global control very fast, as the limit of the zigzag instability 
of solitary band solutions shifts to the left and upwards, and at the same time the tip of the 
cusped region of stable disk solutions retreats to higher values of both fj, and v. 

This picture is in qualitative agreement with the results of the computer experiment 
by Middya and Luss Hl5| , although in the latter work there was no separation between 
characteristic scales of both local variables. In runs without global control, Middya and 
Luss observed, depending on initial conditions, either bands or disks. Although typically 
a large number of either objects were present, they could be seen, essentially, as solitary 
structures, since the diffusivity of both local variables was small, and therefore separated 
kinks could interact only very weakly. With the increasing global control parameter, the 
number of either bands or disks decreased, and some bands disintegrated into a number of 
disks, apparently, due to zigzag instability. 

One should recall that the global control parameter f3\ or /3 2 , as defined above, is pro- 
portional to the number of "solitary" objects N . Therefore decreasing N brings the system 
back into the parametric region where the remaining solitary structures can be sustained. If 
the values of \x and v correspond to the parametric region where both bands and disks are 
stable without global control, bands typically disappear first as the global control becomes 
stronger. Subsequently, the sustainable number of disks decreases, and, eventually, inhomo- 
geneous solutions disappear altogether, in agreement with the observations by Middya and 
Luss. 
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FIGURES 




FIG. 1. Synclinal (a) and anticlinal (b,c) disposition of null isoclines. 
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FIG. 2. The profile of the long-range variable v for a solitary band (a) and disk (b) solutions. 
The value of v is normalized by fx, and the value of the parameter v is chosen in such a way that 
the bounding kink is stationary in the absence of global control. The curves are marked by the 
values of the band half-width or the disk radius a. 



. 005 




1 2 3 4 5 6 

FIG. 3. Construction of stationary solitary disk solutions. Each curve presents the plot of the 



r.h.s. of Eq. (32), and is marked by the appropriate value of the global control parameter fo- 



Solutions are obtained as intersection points with straight lines presenting the l.h.s. of Eq. (32) 
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0.5 1 1.5 2 2.5 3 
FIG. 4. Dependence of the critical value of \i at the threshold of the long-scale concerted zigzag 

instability of a symmetric striped pattern on the half-period L. 
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FIG. 5. Fxistence and stability regions of the solitary band solution in the parametric plane 
(//, u). The solution exists in the region bounded by the diagonal v = \i and the solid line marked 
by an appropriate value of the global control parameter f3\ . Dashed lines denote the limits of the 
zigzag instability at the same values of the parameter (3\. 
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FIG. 6. Existence and stability regions of the solitary disk solution in the parametric plane 
(/[/, v) for the values of the global control parameter = 0, /?2 = 10~ 2 (a) and = 0, /?2 = 10 -3 
(b). Dashed lines denote the instability limits in the quadrupole mode. At = 10~ 3 , the instability 
limit is outside the depicted region. 
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0.8 1 1.2 1.4 

FIG. 7. Comparison of the existence and stability regions of the solitary band (B) and disk (D) 
solutions in the parametric plane (fi, v) for the values of the global control parameter fa = fa = 
(a,b), and fa = 1CT 3 (c). In (a,b), solid lines denote the bounds of existence, and dashed lines, the 
stability limits of both the band and the disk solutions. The region adjacent to the cusp is shown 
on a larger scale in (b) . The solid cusped curve in (c) shows the existence bounds of stable solitary 
disk solutions at fa = 1CT 3 . The dashed lines in (c) denote the existence bounds of stable solitary 
band solutions, and are marked by the appropriate values of fa. 



an 



